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Abstract 

The relativistic Boltzmann equation for a single particle species generally implies a fixed, un- 
changeable equation of state that corresponds to that of an ideal gas. Real- world systems typ- 
ically have more complicated equation of state which cannot be described by the Boltzmann 
equation. The present work derives a 'Boltzmann-like' equation that gives rise to a conserved 
energy-momentum tensor with an arbitrary (but thermodynamically consistent) equation of state. 
Using this, a Lattice Boltzmann scheme for diagonal metric tensors and arbitrary equations of 
state is constructed. The scheme is verified for QCD in the Milne metric by comparing to viscous 
fluid dynamics. 



1 



I. MOTIVATION 



The Boltzmann equation is a tool that has proven to be very useful in many different areas 
of physics. Despite its usefulness, there are some properties of the Boltzmann equation that 
are not optimal for modelling physics systems. In particular, for a single particle species the 
equation of state is fixed by one parameter alone, namely the particle's mass. Since in the 
limit of a small particle mean free path the Boltzmann equation describes fluid dynamics, 
this implies that the equation of state for the fluid hence described is unchangeable and 
(typically) not realistic. This is a problem in particular for the so-called Lattice Boltzmann 
Approach to fluid dynamics jl-3], where the Boltzmann equation serves as a convenient 
algorithm for computing the behavior of fluids. In the non-relativistic context, ways to 
circumvent this problem are known, e.g. a modification of the equilibrium distribution 
function, modifying on.y the pressure components or introducing a new force term ffl, 
but it is not obvious how to generalize those to the relativistic case [9|, |10| . 

On the other hand, it is known that particle masses change when considering a heat bath: 
for instance, photons acquire temperature-dependent masses in a plasma, which leads to a 
corresponding change of the plasma equation of state In high-temperature Quantum- 
Chromodynamics (QCD), these medium-dependent quasiparticles have been successfully 
used to model the QCD equation of state jl2l ]. Is it thus possible to write down a 'Boltzmann- 
like' equation for a single particle species with a medium- dependent mass that can reproduce 
any thermodynamically consistent equation of state? The objective of the present work is 
to give an affirmative answer to this question by means of an explicit construction. 

Note that in the context of quasiparticle and Nambu-Jona-Lasinio models, essentially all 
the relevant parts of the present derivation can be found 13H17|. In this sense, the results 



presented here are not new. However, as far as I can tell, all published results employ 
multiple species of particles, while the results below are for a single species of a 'virtual' 
particle, and therefore probably computationally cheaper. Also, to my knowledge the present 
work is the first to provide a concrete example for an algorithm outside equilibrium with 
an arbitrary equation of state. The article is structured as follows: in Sec. [TTJ I give a 
textbook-style review of the Boltzmann equation in curved spaces. In Sec. 1111} a framework 
for arbitrary non-ideal equations of state is set up and subsequently tested for the case of 
QCD at high temperature. In Sec. IIV} a relativistic lattice Boltzmann scheme for matter 
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with a non-ideal equation of state in curved spacetime is given, with the particular example 
of QCD in a Milne spacetime that may be of relevance for high energy nuclear collisions. 
Finally, I conclude in section |V] 



II. BOLTZMANN EQUATION IN CURVED SPACE: A REVIEW 

This section gives a text-book style review of the Boltzmann equation in curved space, 
introducing the usual particle current and energy-momentum tensor. Expert readers may 
want to skip this section and read on in Sec. IIHI 

The Boltzmann equation specifies the evolution of the single particle distribution function 
/(X 11 , P M ), which is dependent on space-time X M = (t, x) and four-momentum P M = (E, p). 
If collisions are absent, but forces such as gravity are present, particles are assumed to prop- 
agate along geodesies which can be parameterized by an affine parameter T. Accordingly, 
the single particle distribution / does not change along geodesies, 

df dt df dx df dP a df 

Multiplying with the mass m one can recognize mdt/ dT = E, mdx/ dT = p, the energy and 
momentum of a relativistic particle. When re-instating collisions, particles will no longer 
follow geodesies, so df /dT will no longer be vanishing. Hence in the general case one has 

P»dJ + F a dWf = -C[f], (2.1) 

where C[f] is the collision term and F a = tti^^- the force felt by individual particles. For 
gravity, the force is given by F a = —V^ V P^P V where r° are the Christoffel symbols that 
are calculated as derivatives of the underlying metric tensor g^. For electromagnetism, the 
force is given by the Lorentz force F a = qF a ^P^ where F a @ is the electromagnetic field 
strength tensor that can be specified in terms of electric and magnetic fields, and q is the 
particle's charge. 

Including both the gravitational and electromagnetic force terms, let us now take an 
integral moment of Eq. (12. ip with weight 

J d X = J ^ v / ^2e(p°)(27r)5 [g^P v - m 2 ) , (2.2) 



(2nY 

where for clarity d 4 P = n^=o dP 1 * and 9 denotes the Heaviside step-function, g denotes the 
determinant of the metric tensor g^ v and I have adopted the 'mostly-minus' sign convention 



for the metric. The delta-function in d X places particles on the mass shell and the step- 
function picks out positive energy states. Apart from the appearance of p°, which could be 
replaced by a scalar product with a future pointing four-vector, this form of d X is Lorentz 



covariant (cf. |l8j]). Using d^^g = ^/=g~T" and d x g pu = ^ p x 9 pu + V p Xl/ g pix one has 



V M [2Q(p°)5 (P 2 - m 2 )] = 2e(p°)5' (P 2 - m 2 ) 2P a P p T^ 



where V M denotes the (geometric) covariant derivative. Rewriting 2P$ (P — m ) 

-)(p 



dfb (P 2 - m 2 ) and using partial integration one finds 



J d X P»dJ = V,J d x P»f + J d X T^P a P»dff. 

Also, it is straightforward to show that / d X F a/3 'Ppd%' f = via partial integration and the 
fact that F a/3 = —F^ a . Hence the Boltzmann equation implies 

V„/ d x P»f = -J d X C[f]. (2.3) 

Defining the particle number current as = J d X P^f one finds that the Boltzmann equa- 
tion implies the covariant conservation of particle number, V^N^ = 0, if / d X C[f] = 0. 
Taking the integral moment / d X P v of Eq. ( 12. ip , one finds 

V M J d X P'P u f - qF uP J d X Pp = -J d X P u C[f] . (2.4) 

For uncharged particles (q = 0), and defining the energy- momentum tensor as T^ v = 
/ d X P fl P u f, the Boltzmann equation implies covariant conservation of energy and momen- 
tum if 

Jd X P u C[f] = 0. (2.5) 

I will assume the collision term to fulfill Eq. ( 12. 5 p for the remainder of this work. For charged 
particles, the Boltzmann equation implies 

V^ = qF^N P) 

or the change of energy and momentum being caused by the Lorentz force for a current 
J/3 = qNp. For the remainder of this work, I will deal with uncharged particles (q = 0). 
However, the generalization to charged particles should be straightforward. 
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A. Equation of State for Uncharged Boltzmann Gas 



In equilibrium, the energy-momentum tensor is given by ideal hydrodynamics, 

T% = eU^U u - pA"" , (2.6) 

where is the fluid velocity obeying U 2 = 1 and A MV = — U tl U u . The equilibrium 
energy density e and pressure p of the system are related by the equation of state. Since 
Eq. (12. 6p must correspond to the particle's energy-momentum tensor in equilibrium, one has 

e = UpUrTg = Jd x (P^,) 2 / eq , p = = -\jd X [P 2 ~ (P^) 2 } f eq , 

which may be conveniently evaluated by performing a Lorentz boost to the frame where 
P^U^ = p° (recall that d\ is Lorentz covariant). 

Let us now consider a specific equilibrium distribution function for a system of uncharged 



particles (cf. |18|). 

f eq (X a } P«) = Zx exp [- ( PaU ^~ fl )] , (2-7) 

where Z denotes the number of degrees of freedom and [i and T are the chemical potential 
and temperature, respectively. In Eq. (12 .7p . U a is a macroscopic velocity that can be identi- 
fied with the fluid velocity in Eq. (I2.6p . In this case, e,p and the number density n = U^N^ 
may be evaluated as 1 

Ze^ /T m 2 T ( [m\ fm\\ Ze^ /T m 2 T 2 fm\ 
£= 7T-, 3TX 2 - ) + mK x (- , p= — K 2 (-), n = p/T , 



2tt 2 V \TJ 1 \TJJ ' * 2tt 2 \T, 

by using the identity J,™(x 2 - m 2 ) n+1 ' 2 e- x l T = (2n + l)\\K l+n (m/T)(mT) n+l for modi- 
fied Bessel functions K a . It is straightforward to show that these results obey the basic 
thermodynamic relations 

e + p = sT + pn , de = Tds + pdn, (2.8) 

where s denotes the entropy density. From the equation of state, an interesting quantity 
to calculate is the speed of sound squared c 2 = dp/de. For illustration, at p = it can be 
calculated from the above expressions as 

m K2(m/T)\ 1 



c 2 s (T,p = 0) 



T K 3 (m/T) 



1 Note that this definition of n corresponds to dp/d^i. To see this, first go to the local rest frame where 

cq/Op = -<9o P) /cq- 



P 2 — (P a U a ) 2 — p 2 and then rewrite df cq /dp = — 9q / eq . Integrate by parts and rewrite 2p 2 S'(P 2 — m 2 ) 



p l d\ p ^5{P 2 — m 2 ). Another integration by parts then gives dp/dp = U^N^. 
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which increases monotonically with temperature from zero to 1/3. Also, the relation p = nT 
is the equation of state of an ideal gas. Clearly, non-ideal equations of state with a non- 
monotonic behavior of c s or p ^ nT are not describable in this framework. 

In particular, note that changing the behavior of the equilibrium distribution function 
/ eq will not change the relation p = nT, and hence does not provide the freedom needed to 
describe a particular non-ideal equation of state that is dictated by nature. 



III. NON-IDEAL EQUATIONS OF STATE 

As shown in the preceding section, the Boltzmann equation (12. ip for a single uncharged 
particle species leads to equations of state that depend only on one parameter, namely the 
particle's mass. In order to describe arbitrary equations of state with a single uncharged 
particle species I therefore want to consider temperature (and density) dependent masses 
m — > M(T,p), motivated by the fact that in a plasma at high temperature or density this 



approach is physically sound [ll|] . The particles described by the Boltzmann equation should 
then be regarded as virtual or 'quasi'-particles, but for sufficiently non-ideal equations of 
state, they will no longer correspond to any real excitations found in nature. However, 
the virtue of introducing this virtual particles will be that no long-range forces or particle 
mixtures will be necessary to describe the macroscopic system dynamics. 

One immediate problem that arises when considering medium-dependent masses is that 
thermodynamic consistency is no longer guaranteed. Specifically, basic thermodynamic re- 
lations imply that 

fin f)r) 

(3.1) 



€ + P = T &T 



dp 
+ ^ 7T 



T 

which would be violated when inserting m — > M(T, p) in the formulas from Sec. Ill Al To 
fix thermodynamic consistency, I propose the following alternate definition for the energy- 
momentum tensor: 

= J d X P»P u f + B(T, p) gT , (3.2) 

where B(T, p) is a function that will be determined by requiring thermodynamic consistency 
in equilibrium, cf. (13.11) . Calculating energy density and pressure from (13.21) . one finds that 
B(T,p) drops out in e + p and that thermodynamic consistency requires 

= dB + l - J d X f eq dM 2 , (3.3) 




FIG. 1: Left: Results for M(T) when fitting the entropy density from lattice QCD collaborations 

□ □ 
(hotQCD [19( and Wuppertal-Budapest (WB) [20|], respectively) or an interpolation from hadron- 



resonance gas to perturbative QCD (Laine/Schroder, Ref. 



2l\\). Right: quality of the fit (symbols) 



when comparing the trace anomaly e — 3p to the original lattice QCD results (full lines). 

where I used 2p 2 5'(P 2 — M 2 ) = —p l d^S(P 2 — M 2 ) and integration by parts. 
Considering the concrete example (12. 7p . one has explicitly 

e = [SK 2 (f ) + f Kl (f )] + B(T,A , 



P 



Ze „/T M 2 T 2 
2tt 2 



n 



Ze^/ T M 2 T 
2tt 2 



K 2 (f ) . 



(3.4) 



A. Example: QCD at Small Densities 

Let us consider the above construction for the QCD equation of state at zero baryon chem- 
ical potential. In order for the Boltzmann energy-momentum tensor to correctly reproduce 
the high temperature limit of QCD with iV c = 3, Nf = 3, one has to set 



U{Nl - 1) + 7i\yv 



180 

Then, one can determine M(T) by inverting ^ = -^M 3 K 3 (jf) = siqcd, where siqcd 



may be obtained from the lattice QCD results (cf. |l9l. |20|. both Nf = 3) or an interpolation 
from hadron resonance gas results to perturbative QCD (cf. [21]], Nf = 4). Thermodynamic 
consistency requires — Z ^ 2 T Ki ^ = dB }^ , which can be solved for B{T) numerically 
by integrating up from small temperatures where B ~ 0. The resulting fits for the masses 
and the quality of the fit for the quantity e — 3p for three 'physical' QCD equations of state 
are shown in Fig. [TJ 
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B. The 'Boltzmann-like' Equation 



The modified energy-momentum tensor (13. 2p is no longer expected to correspond to a 
moment of the Boltzmann equation (12.11) . because of the extra term in ( 13. 2p . However, 
one can ask if there is a modified 'Boltzmann-like' equation that will give V^T^ = (for 
uncharged particles). Inverting the steps leading to this equation in Sec. [Til an d rewriting 
P^5'{P 2 - M 2 ) = ffi p) 5(P 2 - M 2 ) and integrating by parts I find that 



V M T^ = J d X P u [p^ - r x a() p a pPd? + \d^M 2 dl P) 

where the term involving B(T,[i) cancels if 



f = 



= dB + -[ d X fdM 2 . (3.5) 

Note that this is the same as the thermodynamic consistency condition fl 3 . 3 j) . except that it 
is promoted to hold also out of equilibrium. As a consequence, the 'Boltzmann-like' equation 

P^f - T^P a P^d?f + \d,M 2 dl p) f = -C[f) (3.6) 

is guaranteed to reproduce a conserved energy-momentum tensor that allows arbitrary equa- 
tions of state parameterized by medium-dependent masses M(T,fi). Note that the 'new' 
term ^d jjL M 2 d^f precisely corresponds to the result found when deriving the Boltzmann 



equation from quantum field theory (cf. 22]). Finally, repeating the steps in Sec. |TTJ one 



finds that for J dxC[f] = 0, Eq. (ESJ) leads to 

V„jV" = 0, N" = J dxP^f, (3.7) 
so that the current is formally unchanged when considering medium- dependent masses. 



C. On-shell-ization 



For some applications, it is useful to explicitly perform the dp integral of the Boltzmann- 
like equation. The reason is that if one is interested in moments of the Boltzmann equation 
with respect to the integral measure dx, this will allow one to work with a distribution 
function / that then only depends on p rather than the four momentum P^. When dis- 
cretizing momenta (see below), one thus only has to deal with three dimensions rather than 
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four. Note that any factors of p° that one may have wanted to include before integration 
will simply turn into multiplicative factors of E = <J {.—9ijV % P + M 2 )/ 'goo because of the 
delta-function that is part of d\- Defining 



/(X",p) = / dp 2p°Q(p°)6(P 2 - M 2 )f(X»,p°,p) 



(3.8) 



which is in accordance with 23] up to a factor of goo, one can integrate Eq. ( 13.61) with 
J dp°2Q{p°)5(P 2 - M 2 ), finding 



% [P^\~ r a A P) [ P a P^ ) + ^m 2 d\ p) ± + 2T a P^ = ~C[f) , 



./ 



/ 



/ 



E 



E 



E 



E 



(3.9) 



where P M here is to be understood as on-shell momentum, P^ = (E, p). In terms of /, the 
particle current and energy-momentum tensor are given as 



J (2k¥ V 9 E' J 



f 



(2k) 



d 3 p 
(27p 



- g P»P»L + B(T,n)g>". (3.10) 



As a specific example (cf. 231]) . consider a metric with a line element of the form ds 2 = 
dt 2 - R(t) 2 (dx 2 ). Then r a ^P a P^ = 2Ep i R'/R , T^P^ = 3E R'/R and Eq. fl3U) becomes 

As another example, consider the line element ds 2 = dr 2 — dx 2 — dy 2 — r 2 dY 2 (Milne 
metric). Then one has r P a P» = -2Ep Y g iY T , PP" = E/t, so that one finds 



2p Y E 



P"0,J + \diM 2 dy = -C[f] 



(3.11) 



IV. LATTICE BOLTZMANN-EQUATIONS 

The main idea behind Lattice Boltzmann equations is to have a minimum sampling of 
momentum space given by a discrete set of N vectors P% with n — 0, . . . N — 1 such that the 
conservation equations for the current and energy-momentum tensor are reproduced exactly. 
For maximum efficiency, one uses a linear ansatz for the collision term 

C[f] = ^ (f - Q (4.1) 

with tr the relaxation time. Taking the first two moments of the Boltzmann equation this 
leads to the conservation of the current and energy momentum tensor provided that 

U,T^ = U,T% = e(T, n)U" U,N» = U,N» q = n(T, ») , (4.2) 
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where the equilibrium energy and particle densities are given in Eqns. ( 13. 4p . The function 
B, which is required to always match its equilibrium value, is determined from Eq. (13. 3p . 

Before discretizing momentum space on a lattice, it is instructive to first consider the 
shear and bulk viscosity coefficients that the collision term (14. ip corresponds to. 



A. Chapman-Enskog Expansion 



In the hydrodynamic (close to equilibrium) limit, the particle distribution function can 
be expanded around equilibrium in powers of space-time gradients, 

/ ./;, • /i •/,••••• 

where f\ is of first order in gradients, ji of second order, and so on. In the absence of 
external forces (r^ = 0), the Boltzmann Equation (I3.6P with the collision term ( 14. ip can 
then be solved iteratively in powers of gradients. Specifically, to first order in gradients one 
finds 



fi 



P-U 



'fj,Jeq -r ^li lyJ - ^(pjjeq , 

which can be evaluated easily using / eq = Ze~ P U / T . Since a small-gradient expansion 
corresponds to an expansion around ideal hydrodynamics, we may use the equations of ideal 
hydrodynamics to simplify the above equations. Specifically, for a metric signature of the 
form H one has (c.f. |25| ) 

D\ns = -V - U , DU a = c 2 s V" In s , 

where c s is the speed of sound and D = U^d^ and V Q = A a/5 ^, A a/3 = g a/3 — U a U^. Using 
the thermodynamic relation ^ = c 2 s T = and consistently ignoring higher order gradient 



term corrections, one finds 



fi = f. 



tr 
sq P-U 



rp ®HV 



P 2 - (1 - 3c 2 )(P ■ U) 2 - 3c 2 s MT-^) V • U 



3T V" v ~ "~ s/v ~ ~' ~~ s "~ dT 
where a^ u = V ^U v ) — |A Mi/ V ■ U. Decomposing the full energy momentum tensor into 

= J dxP^P v f + g^B = eU^U" - PAT + tt^ + A^IT , 

where 7r AtI/ = 2r]cr tlu and II = ( V ■ U, one identifies the shear and bulk parts of the dissipative 
tensor with 

^ u = T<^> = J d x P <fl P u> f , n = 1 -a^ + P = \J dx^P^h ■ 
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7T f 



Note that / d\(P • U) 2 f\ = because of Eq. (12. 5p . Using usual the decomposition of the 
integrals in a tensor basis spanned by U tl U u and one finds the shear and bulk viscosity 
coefficients from ir^ u and II as 

V 15T J (2tt)3 P 2 ' 1 } 

t TR r f d3p 2 -M 2 + (l-3c 2 )E 2 + 3c 2 s MTdM/dT _ E/T 

^ 9T J (2vr)3 P E 2 6 ' 1 } 

where E = a/M 2 + p 2 . After a little bit of algebra it is possible to show that in the 



massless limit rj = tr ^"t^ (c.f. [26j), while for constant masses rj = ^ J T c?T" (e + P), 
C = (-3c£r(e + P) + 5 ; T dT' (e + P)) . No simple formulae seem to exist for medium- 
dependent masses. Note also that these results differ from Ref. j^, Q (and many others 
using the Israel- Stewart ansatz) because non-linearities where not taken account there prop- 
erly. 



Pushing the Chapman- Enskog expansion to second order or following Ref. [28| would be 
desirable to extract all the second-order hydrodynamic transport coefficients 25]. While 
this is left for future work, it is possible to extract the value of the hydrodynamic relaxation 
times for the shear sector, t w . Identifying the relaxation time with the coefficient that is 
multiplying —U a d a {j]a^ v ) in 7i^ u , one finds after a little algebra 

from the derivative of f\. Note that it is therefore possible to use the same value of T n 
in numerical simulations using the Lattice Boltzmann with medium- dependent masses and 
second-order hydrodynamics. 



B. Lattice Boltzmann with Medium-Dependent Masses 

In the following, I will present a minimal set of vectors P£ that is usable for a general 
relativistic Boltzmann equation with medium dependent masses, albeit only for metric ten- 
sors that are diagonal. The scheme is constructed by noting that the (on-shell) equilibrium 
distribution function for a Boltzmann gas can be expanded as 

oo 

/ eq pP, p) = e^~ E ^ T £ (p^/T) n jn\ , 

n=0 
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and I recall the definition of E given in Sec. IIII CI If the metric is diagonal, one may 
rescale the space-like momentum components such that —gijP 1 ^ — >■ bijfi 1 !? = |p| 2 - (Note 
that this rescaling also changes the form of the Boltzmann equation.) Now p/|p| is a unit 
vector that may be parameterized by spherical coordinates (angles 0, 9). Setting furthermore 
|p| = Msinh£ (implying E = M cosh £ / y^oo) one has the parametrization 



P M = M sinh £ (cotanh£ / y/goo , sin 9 cos / y/—gn, sin 9 sin 0/ V — ^22, cos / V — fe) 
for the momentum in terms of the variables £, 0, 0. Therefore one has 

n=0 

where v 1 = p*/M = (sinh £ sin cos 0, sinh £ sin 9 sin 0, sinh £ cos 9) and o' n '(X' 1 ) are some 
coefficients that are only space-time (but not momentum-) dependent. Powers of the veloci- 
ties 2 v may be be represented using the polynomials P^, in (v/|v|) that are orthogonal with 
respect to the angular integral dtt (see Ref. [lo| for details). 

This then motivates the ansatz for the general distribution function: 

00 00 / \ 

/(I^^.^e^^^E^l 71 sinh"^ (fc) (coshOaSi(^) ? (4-5) 

n=0fc=0 \I V I/ 

where To, Mq are some reference temperature and mass, respectively, and R k are polynomials 
of degree k that will be defined below. In practice, the infinite sums above are truncated at 
some finite order. Furthermore, it turns out that for any even n, the sinh£ terms may be 
represented by the sum over polynomials R{k)i so Eq. (14. 5 j) may be modified such that there 
is a single inverse power of sinh £ for every n odd. Replacing continuum momenta P M by a 
discrete set requires the condition that the integrals in Eqns. (I3.10p are represented exactly. 
For the angular integrals, this requirement is identical to that of massless particles discussed 



in Ref. [10[. Note that Eqns. (I3.10p are then evaluated for fixed values of £ rather than 
fixed E, p, implying another change in the Boltzmann equation coming from the space-time 
dependent masses. For convenience, a concrete example will be given below. 



Note that v really is equal to velocity times the Lorentz factor. 
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C. Deriving the Momentum Lattice 



Let us quickly review the derivation of the discrete set of momenta: the angles can be 
found from the requirement that 

c2tt 



rZTT -77- 

/ #(sin0) a (cos0)^- a - 1 = — E(sin0O°(cos0i) ' ' 

JO A rh i 



where a is assumed to be a non- negative integer smaller than 2N^ — 1. Namely, the above 
integrand can be recast as a Fourier series involving as highest harmonics cos[(AT^ — 1)0] 
and sin^iV^ — 1)0]. Exact representation of the integral as a sum is possible if the angles 
are chosen as the nodes of functions orthogonal to the integrand. For the Fourier series 
above, there are actually two sets of orthogonal functions: cos[A^,0] and sin[i\^,0]. Choosing 
sm[N^(p], the nodes are given by = 0; = j^-, I = 0, 1, . . . 2N ( j ) — 1, which fixes the set for 0. 
For the discrete set of angles 9, note that the integrands Eq. (13.101) only depend on 9 through 
P M , so that sin 9 always comes with either cos or sin 0. Since any odd power of cos 0, sin 
integrates to zero, any non-vanishing contribution must involve sin 2 9 = 1 — cos 2 9. Hence, it 
is sufficient to consider only integrands with powers of cos 9, which may be recast as a sum: 



f d(cos 9) cos 2 ^- 1 = $>J cos 2 ^- 1 9j , 

1 j 



where the discrete angles 9j are given as the roots of the Legendre polynomial 
L^ e (cos 9j) = 0, j = 0, 1, ...Nq — 1 and the weight factors W® are given as Wj = 
2/ [(l-cos 2 %) (L^(cos^)) 2 . 

Similarly, the integrands in Eq. (13.101) then only depend on cosh£ and sin 2 £, since any 
odd power of sinh£ would have integrated to zero already. Therefore, it is sufficient to 
consider integrals of the form 

/ dee-*^«sinh 2 ecosh^- 1 e = E4(^)cosh 2;v e- 1 (e fc («b)) , z = M /T , 

J ° k 

and the nodes and weights w^. are calculated from the set of polynomials which are 
orthogonal on J G?£e _z ° cosh £ sinh 2 £. Specifically, one finds 

Ro(£) = l, i?i(0 = coshe-fg, 
K 2 {l;) - cosn £ + , o(z 3^3 +8Zo ^ o ^ 1+1 4 ZoX() ^ +2(2 _^ ) ^3_ 2 3^^i(U Z()Kl , etc. 
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Hence the discrete values are calculated from (^k) = and the weig hts w[ fulfill 

A/f — 1 / \ A/f-1 

E = ^^S E vlRmitk) = 0, ro=l,...JV € -l. 

fc=0 Z ° fc=0 



Thus, one finds the following representation of the momentum integrals: 

-JO roo A 7V 5 -liV fl -l2^-l 

/ — / desinh 2 ^/(e,0»= E E E w kj f{t k ,e j ,<j>i) = Y,wnf(Pi:), 

J 47r Jo fc=o i=o /=o n 

with the weights w k j = e z ° cosh ^ k wlwj /(AN^) = io n , and where the collective index n runs 
over all discrete momenta P£ constructed from the ensemble 4>i,9j,^k- 

D. Lattice Boltzmann Algorithm for Milne Spacetime 

In this subsection I give a detailed construction of a lattice Boltzmann algorithm with non- 
ideal QCD equation of state in an expanding spacetime with ds 2 = dr 2 — dx 2 — dy 2 — r 2 dY 2 
(Milne). For simplicity, I will limit myself to neglecting space dependencies, which are 



algorithmically easy to program (cf. [10| for a practical example). 

Starting with the Boltzmann-equation (13. 1 If) for the on-shell distribution function, let us 
first rescale momenta p Y = p Y /r so that E = 5ijp l fP . Next, replacing p l = M(r)v* and 
using Eq. (14. ip . Eq. (13.1 ip becomes 



dj - d T In M p • %)/ = J—h* . (4.6) 



cq 

T - - - T R 



Neglecting space dependencies, one has cylindrical symmetry and hence / = /(r, £,#). 
Therefore the ansatz for the general distribution function can be simpler than (14.51) . namely 

oo oo 

fiX*, e, 6, 0) = e -^ocosh<; J- E L n (cos 0) %)(coshO a^ k \r) . (4.7) 

n=0 fc=0 

Using the nodes and weights from Sec. IIV C\ one immediately finds 

« (m ° = J^J ^sinhe 2 L m (cos^)^(coshO/ = ^y^E^ L^/| p , , 

where 1(1) = J d£ sinh 2 £i? 2 (cosh£). 

I will not consider conserved particle number, so the only quantity of interest is 

M 2 

T " V = V - 2 Y.^nf{P^P u n +B{T)g^. (4.8) 
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Since all spatial dependencies have been neglected, the fluid velocity is trivial, = (1,0), 
and the equilibrium energy density is given by the 00 component of Eq. (14. 8 p . 

A change in the energy density can be calculated via the Boltzmann equation ( 14. 6ft : 

/ /eq 



d T e = — Vw n cosh^ 



p -d?f 

T 



where I used partial integration and the identity (13. 5p . which becomes 

<9 T M 4 



J2w n f(PX) + d T B(T)=0. 



(4.9) 



(4.10) 



More specifically, a lattice Boltzmann algorithm may hence be constructed as follows: a 
valid initial condition at time r consists of specifying / = / cur and an initial temperature 
and particle mass, T cur , M cur . Then, make a prediction of the change in distribution function 
and (logarithm of) mass: 



pred 



'5 1 ! | SlnM pTed ~ _ fen 

t St ~ t r (T c 



feq (-^ci 



<)lnM pred = - — 



J2w n \S x /t- 



/cur /eq (^cur ) 
Tr(Tcut) 



St 



St 



(4.11) 



T cm . S w nf cur 

where S\ and $ 2 are representations of the momentum derivatives p Y and p • d^f in the 
form of (14. 7\\ with coefficients 
2m + 1 



ml 
b 2 



1(1) 

2m + 1 



E w ™/c U r [(P m R[ cosh f - ffc) cos 2 tanh £ 2 + #z (> m + P^(cos - cos 3 

1 



PmRl 2 + 



cosh 2 £ 



+ P m -R; cosh £ tanh £ 



respectively. Via Eq. (14. 9p . this leads to a prediction for the new temperature T pre d. These 
predictions are then corrected using the trapezoid integration formula 



Sfo 



St 



drf +0 T f 



-0(St) 3 , SliiM a 



St 



d T \n M\ t + d T \n M\ t+St )+0 (St 



where the values at time r + St are calculated using Sf prec \ and SliiM prec \. Note that the re- 
sulting mass M corr does not necessarily correspond to the equilibrium particle mass M(T new ). 
As a consequence, I use / eq = Ze~ V M ^ +M ^°" sinh2 £/ T for the equilibrium distribution func- 
tion in the algorithm. The above steps may be repeated to solve the Boltzmann equation 
( 14. 6p for arbitrary times. The resulting algorithm leads to time integrated quantities that 
are accurate to 0(St) 2 (cf. |29|). 
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E. Results for Milne Spacetime 



In this section I provide tests of the above Lattice Boltzmann algorithm by comparing 
results to viscous fluid dynamics for the QCD equation of state of Ref. 2l| and a Milne 
metric. The fluid dynamics equations for the energy density and quantity $ = Ty — p fulfill 
the coupled equations [24] 



e _i_ p $ $ 4r? 4$ Ai 

d T e = - e -^- + -, d T $ = -- + -^-^--V> 2 , (4.12) 



where tv is the relaxation time and Ai is a self-coupling parameter. While I found r n = tr 
in Sec. IIV A\ Ai is currently not known. For simplicity, for the hydrodynamic calculation I 
will use the values r n = and Ai = ^r]r n that are reported for the massless gas case [26 ]. 

One should keep in mind that — since the correct values for r n , X± will differ from this 
choice in view of the findings in section IIV AI — this implies that the hydrodynamic and 
Lattice Boltzmann results will not agree in practice. Note, however, that there is another 
issue that prevents perfect agreement between (second-order) hydrodynamics and Lattice 
Boltzmann theory even in principle: the reason is that, even if one were to use the same 
second-order transport coefficients in a numerical second-order hydrodynamics and a Lattice 
Boltzmann solver, the two would still disagree because of the different third order gradient 
terms. However, for all practical purposes when hydrodynamics itself can be considered 
applicable, the difference between the two numerical schemes could be considered small. 

For the Lattice-Boltzmann framework, the QCD equation of state is parameterized as in 
Sec. IIII AI with a reference value Mq/Tq = 1 for (14. 5p . Note that this reference value corre- 
sponds to a reference temperature of T ~ 0.82 GeV, This means that results will be most 

n 

accurate for this temperature (cf. [10]), and in particular will break down if applied to prob- 
lems involving fluid cells with temperatures exceeding two times this reference temperature. 
The results at temperatures different than this reference value can be improved by increasing 
the value of N^, but in practice I find that Ag > 3 gives adequate results. In Fig. [21 I show 
the temperature evolution in viscous fluid dynamics and the above lattice Boltzmann algo- 
rithm for N% — 5, Nq — 5 and 5r = 0.2 fm/c. As can be seen from this figure, the non-ideal 
equation of state time evolution in fluid dynamics for rj/s = 0.5 is described reasonably 
well throughout the whole simulation time, even though it differs markedly from the ideal 
equation of state time evolution (shown in Fig. [2] for comparison). Overall, the algorithm 
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Temperature 



Scaled Temperature 




T [fm/cj i [fm/c] 

FIG. 2: Temperature evolution for various viscosities, from fluid dynamics ('VH') and Lattice 
Boltzmann ('LB') respectively. Results are for a QCD equation of state except for the ideal equa- 
tion of state evolution (left plot). Transport coefficients for VH are only approximate, so perfect 
agreement is not expected. Right: results normalized with respect to the ideal fluid dynamics to 
highlight differences. The r]/s = oo results are obtained by setting tr = oo (free-streaming) and 
updating the medium-dependent mass according to the change in energy density. These results do 
not correspond to an actual physics situation and are presented for illustrative purposes only. 

seems to perform rather well and provides a concrete example for simulating the Boltzmann 
equation for a system with non-ideal equation of state and non-Minkowski geometry. 



V. CONCLUSIONS 



In this note I have set up a general relativistic transport equation for a single species 
of uncharged particles with a medium-dependent mass. This 'Boltzmann-like' equation 
allows for a conserved particle current and energy-momentum tensor. The latter explicitly 
allows arbitrary (thermodynamically consistent) equations of state when using the medium- 
dependent mass as fitting parameter. I expect this formulation to be useful for relativistic 
fluid dynamics simulation in the Lattice Boltzmann framework with arbitrary equations of 
state. Possibly, it can also have relevance for non-relativistic computational non-ideal fluid 



dynamics (cf. 30] ). Furthermore, it can be applicab 



dynamics with non-ideal equations of state, cf. 



31 



e in the context of simulating parton 
32l | or quasiparticle models of QCD [16 1. 
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